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The collective ability of organisms to move coherently in space and time is ubiqui- 
tous in any group of autonomous agents that can move and sense each other and the 
environment. Here we investigate the origin of collective motion and its loss using 
macroscopic self-propelled Bristle-Bots, simple automata made from a toothbrush 
and powered by an onboard cell phone vibrator- motor, that can sense each other 
through shape-dependent local interactions, and can also sense the environment 
non-locally via the effects of confinement and substrate topography. We show that 
when Bristle-Bots are confined to a limited arena with a soft boundary, increas- 
ing the density drives a transition from a disordered and uncoordinated motion to 
organized collective motion either as a swirling cluster or a collective dynamical 
stasis. This transition is regulated by a single parameter, the relative magnitude of 
spinning and walking in a single automaton. We explain this using quantitative ex- 
periments and simulations that emphasize the role of the agent shape, environment, 
and confinement via boundaries. Our study shows how the behavioral repertoire of 
these physically interacting automatons controlled by one parameter translates into 
the mechanical intelligence of swarms. 
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Collective behavior is ubiquitous among living organisms: it occurs in sub- 
cellular systems, bacteria, insects, fish, birds and in general in nearly any group 
of individuals endowed with the ability to move and sense (Miller 2010, Vicsek & 
Zafeiris 2012). Recent studies of collective behavior have focused on the mechanism 
that triggers the switch from disordered to organized motion in a swarm (Vicsek 
& Zafeiris 2012, Vicsek et al 1995, Gregoire & Chate 2004, Ballerini et al 2008, 
Leonard et al. 2012, Buhl et al. 2006), and its implications for artificially engineer- 
ing these strategies in robotic systems (Mallouk & Sen 2009, Rubenstein et al. 2011, 
Mellinger et al. 2010). For example, in social insects, such as the agarophilic desert 
locusts, the transformation from solitary to social behavior arises as a consequence 
of proximal tactile interactions that are density controlled (Buhl et al. 2006). Exper- 
iments on the claustrophilic termites, Macrotermes michaelseni which are used to 
living in confined spaces, have demonstrated the existence of a variety of collective 
behaviors such as coordinated circulation and arrest or stasis in a closed confined 
geometry. These different behaviors may be triggered by varying the density of 
the colony and disturbing it through external stimuli (Turner 2011). Understand- 
ing how these biological behaviors arise from a mechanistic perspective has been 
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Figure 1. (a) A collection of the BBots used in the experiment, (b) Schematic of an 
individual BBot. A plastic chassis is connected to a pair of toothbrushes via a slanted 
wedge. An eccentric motor is positioned on the top side of the device and is powered by 
a VARTA rechargeable button-cell battery. 



difficult given our primitive experimental abilities to probe the neuro-ethology of 
these complex creatures. Theoretical attempts to understand these behaviors use 
putative models of interactions between organisms as a function of their density in 
periodic domains (Vicsek & Zafeiris 2012), while a practical approach circumvents 
the question of mechanism and implements workable strategies to actively direct 
the collective dynamics of ensembles of agents (Rubenstein et al. 2011, Mellinger 
et al 2010) using feedback control in individual agents (Braitenberg 1984). These 
approaches clarify the common bases at the heart of all swarming behaviors: the 
ability of an agent to move, the ability to sense others and the environment, and 
the ability to respond to both of these kinds of stimuli. 

Here, we probe the transition from random swarming to collective motion and its 
loss using a minimal system composed of self-propelled automatons that can sense 
each other mechanically through contact and interact both with an environment 
of varying topography and with boundaries. Our setting is macroscopic, control- 
lable and especially suitable to investigate the role of the environment in selecting 
and tuning the collective behavior of the group. Unlike experiments on vibrated 
particles (e.g. Narayan et al. 2007, Deseigne et al. 2010), where all particles are 
simultaneously driven using the same source, our agents are autonomous and self- 
propelled, with velocities that are independent, and yet show collective behavior 
even in a small group of individuals in the presence of confinement. 



1. Motion of an individual BBot 

Our experiments were carried out using a custom-fabricated swarm of Bristle-Bots f 
(BBots) (Murphy 2009, Bobadilla et al. 2011), simple self-propelled automata with 
similarities both to natural mechanical ratchets (Kulic et al. 2009) and their arti- 
ficial analogs (Mahadevan et al. 2004). Our system has three controllable features: 

t See also: http://www.evilmadscientist.com/article.php/bristlebot 
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Figure 2. Principle of motion of a single BBot. The bristles act as legs that are periodically 
flexed under the action of the eccentric motor. Over a cycle, the bending and unbending 
of the tilted bristles causes them to slip, resulting in forward motion. 

(i) a tunable ratio of linear speed and rate of turning for individual agents; (ii) 
a collective ability to exert aligning forces and torques on each other by means 
of shape dependent contact interactions; (iii) confinement induced by soft or hard 
boundaries. The design of our BBots (Fig. 1) is optimized to be small, light, stable 
and modular. An elliptical plastic chassis (major axis 7.92 cm, minor axis 1.85 cm) 
serves as a container for a 1.2V VARTA rechargeable battery which can slide inside 
the chassis to adjust the position of the center of mass and thus change the relative 
ratio of translational and rotational speed. The battery is connected to a motor 
(commonly used in cell phones) housed on the top side of the chassis, with a mass 
of 0.5 g and an eccentricity of 0.8 mm, designed to rotate at 150 rounds per second. 
Two rows of nylon bristles, obtained from a commercial toothbrush, form the legs 
of our BBots. The bristles are cut to 5 mm length to prevent tipping without com- 
promising their flexibility and are attached to the chassis via a removable wedge. 
This allows us to control the inclination of the bristles relative to the chassis. The 
total mass of the object is 15.5 g. 

BBots move when the eccentrically loaded motor drives the legs of the machine, 
the bristles, which flex periodically. The bending of the tilted bristles on the sub- 
strate causes them to move more easily in the forward direction relative to the 
rear, leading to a rectification of the periodic driving and thus directed movement. 
Over each period of rotation of the eccentric motor the sequence shown in Fig. 2 
is followed (see also Appendix A and Supplementary Movie 1): 1) the bristles are 
loaded by a force F = Mg + mru'^ {Mg weight of the BBot, m eccentric mass, 
r lever arm, oo angular frequency of the motor); 2) as the eccentric mass rotates, 
the load on the bristles decreases, causing the bristles to recoil; 3) the bristles slip 
forward on the underlying substrate, producing a net displacement of the object. 
To quantify the motion of an individual BBot, we analyzed the shape of a row of 
bristles treated as a single elastic beam subject to a periodic tip load as well as a 
frictional force in the horizontal direction (see Appendix A) and showed how the 
linear velocity of a BBot and its turning rate depend on the design parameters of 
the system. 

In Fig. 3, we show how the motor speed, bristle position, length and angle and 
the system mass leads to changes in the speed of an individual BBot confined to 
a narrow channel to prevent lateral drift. We see that the bristle inclination and 
length have a strong effect on BBot locomotion; increasing the angle a of the bristle 
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with respect to the vertical direction causes the BBots to slow down substantially 
when a is varied from 5° to 30°. The length of the bristles affects the motion of a 
BBot in two ways: longer bristles cause the center of mass to be displaced further 
in each step, leading to a linear increase in the velocity (Fig. lc-3), while short stiff 
bristles lead to a noisier dynamics associated with rebounds and jumps driven by 
the eccentric forcing. Furthermore, because long bristles cause the BBots to spend 
a longer time in contact with the substrate (where the transverse component of 
the eccentric force is balanced by friction), they move primarily along a straight 
line, while BBots equipped with short bristles are prone to move in a circle. This 
sensitive dependence on the bristle parameters allows us to tune the locomotion of 
individual subunits and study its role on the collective behavior of the community. 
In particular, by choosing 5 mm bristles and varying a from 0° (upright bristles) 
to 20°, we obtained two distinct types of individuals: 1) spinners, which are BBots 
with a = 0° and 5° that tend to spin clockwise with an angular velocity of up to 
30 rad/s while moving slowly; 2) walkers, which are BBots with a > 10° that move 
in a straight or weakly curved orbit (Supplementary Movie 2). 



2. The effect of boundaries and topography 

In most prior studies of collective behavior, the boundary is not considered; indeed 
theoretical studies routinely treat only the case with periodic boundaries, while the 
few experimental studies that exist aim to minimize the role of boundaries. In our 
study, the boundaries play a most important role as we now discuss. Our arena 
consists of a circular plate of 44 cm in diameter, with a single BBot taking up 
approximatively 0.8% of the total available area. 
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Figure 4. Experiments showing the interaction of BBot with a soft and hard boundary. 
(a-1) A spinner with a — 5° and (a-2) a walker with a — 10° in a soft-boundary arena. 
(b-1) A spinner with a = 5° and (b-2) a walker with a = 10° in a hard-boundary arena. 
For a soft boundary, a consequence of the shallow bowl-shaped curvature is that the BBot 
is reflected toward the interior. For a hard boundary, the BBot gets aligned with the edge 
of the arena and moves along it. 

(a) Soft boundaries 

We first consider the interaction of a BBot with the boundary that causes it 
to be reflected back from the edge, into the middle, using an arena with a gentle 
upward sloping edge, fabricated by oven- forming an acrylic disc over a frisbee- 
shaped aluminum mold. Here, we see that surface topography plays a role normally 
reserved for the boundary by influencing the motion of a BBot via environmental 
changes. With this soft boundary setup, our BBots either turn back into the middle 
(behavior typical of spinners, i.e. a — 0°, 5°) or they oscillate back and forth in 
a periodic motion that causes them to remain in the neighborhood of a particular 
location at the boundary (behavior typical of walkers, i.e. a = 10°, 15°) (Fig. 4a). 
This pendulum-like effect follows from the fact that the walker's path is never 
perfectly radial, so that as a BBot climbs the edge it also turns sideways. On the 
steepening gradient near the edge, the BBot typically slips backwards, as it rotates 
by about 30°, and Sisyphus-like, tries to climb up the edge again only to be kicked 
back to where it started. These oscillations may be repeated a few times for an 
individual Bbot before it eventually moves back into the center of the arena, and 
then onto another part of the edge where the same phenomena is repeated. Strong 
walkers with a = 20° do not experience the oscillatory motion at the boundary 
because their forward propulsion dominates the role of sideways spinning motion 
and tends to align the BBots to be normal to the edge independent of how they 
initially approach the boundary. 
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(b) Hard boundaries 



To see what happens when we change the environment in which the Bbots 
operate, we replaced the boundary of the arena with a gently curved edge with a 
flat circular disc of the same 44 cm diameter, bounded by a thick (vertical) strip 
of acetate ~ 4 cm high that is held firmly in place by a ring of thick translucent 
tubing. The most salient feature of this hard boundary system is that the boundaries 
are not reflective, so that a BBot that hits the edge will begin to circulate in a 
particular direction around the arena, traveling always parallel to the edge (Fig. 
4b). We observe stable motion in both a clockwise and counterclockwise direction, 
the determining factor being the angle of initial contact with the wall. 

3. Collective behavior of BBots - Experiments 

To study the collective dynamics of the BBots, we use a transparent plate that is 
backlit with a set of neon lamps and allows us to track the BBots with a digital 
camera at 40 fps. The resulting movies were processed with tracking software to 
compute the position, orientation and the translational and rotational velocity of 
each BBot and thus quantify their individual behavior. This allows us to calculate 
the following two order parameters to characterize the collective behavior of the 
putative swarm: 



where TV is the total number of BBots and Vi{t) the velocity of the i—th BBot at 
time t. We see that Vi is the average velocity of the BBots, while V2 is their average 
speed. When they move in a disordered fashion, ^'i ^ (becoming exact in the 
infinite particle limit) and V2 > 0; BBots moving coherently in space have both 
vi ^ and V2 ^ 0, while if a cluster of BBots is dynamically arrested, ^2 ^ 0. 

In our experiments we used a paperboard template that initially arrested the 
motion of the BBots. When this was removed, the BBots moved and eventually 
reached a statistical steady state (Supplementary Movie 3). Alternatively, we pro- 
gressively increased the BBot population from 2 to 16, adding a new one every 30 
seconds and then removing them one at a time, to measure the hysteresis in the 
transition between states of collective behavior. 

In Fig. 5a, b, we show the results of our experiments on the collective motion of 
spinners {N = 7, 24 with a = b°) moving in an arena with a soft boundary for 3 
minutes. Spinners spin rapidly and collide frequently and strongly with each other 
(Supplementary Movie 3); when < 10 (corresponding to 8% area coverage) their 
motion is disordered, with vi < V2 as seen in Fig. 5a-2 and their center of mass 
moves aperiodically as shown in Fig. 5a-3. When > 10 the spinners aggregate at 
the edge of the arena while aligning themselves at an angle to the boundary, and 
start swirling collectively clockwise (the direction of spinning for individual BBots) 
coherently along the edge, as show in Fig. 2b-2. In this case, the order parameter 
vi increases and saturates once the swirling cluster is formed (Fig. 5b-2). 

Walkers at low density have a different behavior than spinners; they move to 
the edge, stay for a while before they turn around randomly, eventually reaching an 
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Figure 5. Experimental realization of collective behavior of Bbots. The columns summa- 
rize the three behaviors observed in the experiments with BBots: (a) disordered (random) 
motion of spinners at low density; (b) swirling motion of spinners at high density, and (c) 
stasis of walkers at high density. (al-a3) Experiments in the random phase; (al) instanta- 
neous position of 7 spinners with a = 5°, (a2) the mean velocity vi and the mean speed 
V2 of the Bbots, (a3) the trajectory of the center of mass of the Bbots in physical space, 
showing random motion. (bl-b3) Experiments in the swirling phase; (bl) Instantaneous 
position of 24 spinners with a = 5°, (b2) vi and V2 showing a non-zero value, (b3) the 
trajectory of the center of mass of the Bbots in physical space, showing the signature 
of the coordinated swirling. (cl-c3) Experiments in the stasis phase; (cl) Instantaneous 
position of 15 walkers with a = 10°, (c2) vi and V2 of the Bbots, (c3) the trajectory of 
the center of mass of the Bbots in physical space, showing no motion, i.e. stasis. 

approximately antipodal point where this behavior is repeated. As the number of 
walkers is increased, they form ephemeral clusters along the edge (Supplementary 
Movie 4) that eventually break up. However, when > 8, clusters of BBots oriented 
perpendicular to the edge form and remain stable, as shown in Fig. 5c {N = 15 with 
a = 10°). This corresponds to the order parameters ^ ^2 ^ 0, and the center 
of mass is essentially stationary (Fig. 5c-2 and 5c-3). In Fig. 6, we show a phase 
diagram that summarizes the collective behavior of BBots confined to an arena with 
a soft boundary, showing disordered motion, swirling and stasis, and highlights the 
hysteretic nature of the transitions between states. For example, once a swirling 
cluster of spinners has formed, it remains stable even when BBots are withdrawn 
from the cluster until N < 6. 
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Figure 6. Phase diagram and hysteresis in collective behavior, (a) A phase diagram sum- 
marizes the three types of collective behavior observed as a function of the number of 
BBots (which is equivalent to their density) and the bristle angle (which is proportional 
to the inverse of the velocity) confined to an arena with a soft boundary, (b) Hysteresis 
diagram obtained by progressively increasing and decreasing the number of spinners in the 
range 2 < N < 16. In the forward portion of the curve (blue) the population transitions 
from disordered to swirling motion for N > 10. After the onset of collective motion, newly 
added BBots are eventually collected by the swirling cluster. When BBots are withdrawn 
from the trailing edge of the swirling cluster the behavior switches from coordinated to 
disordered only when the population is below 6 BBots. 

To understand how confinement and topography lead to these behaviors, we 
first use BBots with an intermediate bristle inclination to observe the assembly and 
disassembly of clusters at a soft boundary; for example, BBots with a = 10° resist 
the sideways motion necessary for swirling but they do not get trapped at the edge 
as easily as BBots with larger a. The result is that clusters can form orthogonal to 
the circular boundary, but if the cluster is too small in number it will eventually 
disassemble due to the growth of coordinated oscillations of the entire cluster (Fig. 
7a). Indeed clusters of four or five bots remain stable for over a minute before 
disassembling. The stasis or jamming region in our phase diagram describes the 
formation of clusters at even higher densities when they become stable over very 
long times. 

In contrast, collective behavior in the presence of a hard-boundary leads to con- 
tact with the vertical wall and aligns the BBots along the boundary thus limiting 
their motion and reducing the interactions between BBots (Fig. 4b). BBots sliding 
along the boundary of the arena eventually form groups due to the small varia- 
tions in the velocity of individuals. However, as the number of Bbots in a group 
increases, it becomes less stable and can abruptly self- arrest. These arrested states 
can take the shape of a half-aster as shown in Fig. 7b when the arena is bounded 
by a vertical, rigid boundary, in contrast with the orthogonally oriented jammed 
structures formed in the soft boundary arena. 

4. Collective behavior of BBots - Theory 

The nature of the collective motion and stasis in our system of confined agents relies 
on the ability of the BBots to march in the direction of their major axis, and rotate 
and align with each other and with the boundary. In order to understand these 
effects quantitatively, we use simulations of self-propelled particles (SPP) consisting 
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Figure 7. An experimental example of cluster assembly /disassembly in the presence of soft 
and rigid boundary, (a) Five walkers with a = 10° initially gather at the soft boundary. 
The cluster, however, starts to oscillate and eventually disassembles, (b) The same walkers 
in an arena with a hard boundary corresponding to a vertical wall jam to form a half-aster 
pattern which nucleates and grows in size until all the BBots in the system have been 
collected by the jammed cluster. 



of two-dimensional ellipses whose center of mass position and orientation Oi are 
governed by the following dynamical system: 

, Ni 

dri , ^ / . X 

= voTii + fei 2^ Fij (4.1a) 

^=c^ + + ^2^M,, (4.1b) 

The first equation describes the over-damped motion of individual ellipses with 
velocity vq along their major axis rii = (cos6>^, sin^^), where Fij is the repulsive 
elastic force between the i—th ellipse and its Ni neighbors, these being defined 
as the set of all ellipses that overlap with the i—th. This force between the i-th 
particles and its Ni neighbors is given by 

F,j=£N,j, (4.2) 

with i a virtual spring lengthy which for the ellipses is calculated from the intersec- 
tions between the two overlapping ellipses as illustrated in Fig. 8. 

The second equation implies that the major axis of each ellipse rotates coun- 
terclockwise with frequency u and can align with its neighbors as a consequence of 
the physical torque due to the contact with the neighbors and is given by: 

Mij = {dij xFij)-z, (4.3) 

where dij is the lever arm of the force Fij exerted by the j — th neighbor on the i — th 
ellipse and z is the unit vector in the z-direction. The constant k2 measures the 
strength of this aligning interaction, while d is a delta-correlated random variable in 
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Figure 8. Schematic representation of the forces between overlapping eUipses. The forces 
are appUed along the direction Nij perpendicular to the line passing through the inter- 
section points A and B of the two particles, at the mid point C. The magnitude of the 
force is controlled by the spring-length i obtained by intersecting the line Nij with the 
perimeter of the region where the two ellipses overlap (shaded in the figure). The overlap 
between the particles is exaggerated in the figure; in the simulations it is very small, so 
that the direction Nij approximates the common normal direction of two convex objects 
touching at one point. 

the interval [— C] represents the noise associated with all the non-deterministic 
factors that affect our system. 

The ellipses are confined to a circular arena of radius R and subject to a non- 
local exponentially decaying torque exerted by the boundary that reorients them 
toward the interior, which reflects the torque produced by the curvature of the 
experimental arena along the edge. The interaction between the particles and the 
boundary takes place through a virtual linear spring acting at the center of mass: 

^boundary ^ _^^(|^.| _ . . ^ ^ ^ ^ 

and a long range torque of the form: 

^boundary _ sin[arctan(y,/x,) - Oi] exp (^lA_^^ . (4.5) 

where ^ is a constant length that can be used to tune the range of the interaction. 
When a particle is in proximity of the boundary [i.e. {\ri\—R)/£^ ~ 1], this boundary 
torque has the effect of rotating the particle toward the interior. The non-locality of 
the boundary torque M^^^^^^^y in (4.5) mimics the distributed gravitational torque 
produced by the curvature of the dish in which the particles move and appears to 
have a significant role for the clustering of the particles at the boundary. Compared 
with its local- analog (i.e. a torque of the same form that acts only when \ri\ > R)^ 
the torque (4.5) has the effect of producing more densely packed clusters of par- 
ticles along the boundary and thus fundamentally changes the collective behavior 
of the Bbots. Being curved, the boundary of the arena has the effect of getting 
the particles to form densely packed clusters. Extending the range of the particle- 
boundary interaction is equivalent to increasing the curvature and thus accentuates 
the focusing effect. This is also the simplest situation where we see how the physical 
environment can control the behavior of these autonomous agents. 

Our system differs fundamentally from those studied in the past by account- 
ing correctly for orientational effects using torque balance rather than an ad-hoc 
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alignment term, while also exploring the role of non-local interactions using topog- 
raphy and finite size boundaries. The dynamical system (4.1) is characterized by 
four dimensionless parameters: the scaled density (j) = /ah (where a and h are the 
minor and major semi-axes of the ellipses), the spinning to walking ratio uja/v^^ the 
orienting parameter /c2a^//ci and the scaled noise parameter (a/vQ. In our experi- 
ments, the relevant experimental variables are the scaled density and the spinning 
ratio, since the orienting parameter and the noise are intrinsic to the shape of the 
agents and the motor characteristic. 

Varying the two relevant parameters in Eq. (4.1) and integrating them numer- 
ically leads to a variety of collective behaviors consistent with our observations as 
shown in Fig. 9 (Appendix B and Supplementary Movie 5). We see that individual 
self-propelled walkers or spinners tend to migrate toward the boundary of the arena 
where they experience a torque that reorients the individual toward the interior. 
At low densities, the primary interactions of these automata are with the boundary 
and so one sees random uncoordinated movements (Fig. 9a). At higher densities, as 
more ellipses simultaneously cluster in the same region of the boundary, the aligning 
force exerted on the ellipses by each other can overcome the action of the boundary 
provided the cluster is large enough. Thus walkers for whom ua/vo ^ tend to 
aggregate into a static structure at the boundary at high enough density (Fig. 9c). 
However, spinners for whom uja/vo > form clusters at the boundary that are 
tilted, and this broken symmetry together with the effect of the weak topography 
(boundary curvature) keeps them confined to the neighborhood of the boundary 
and causes the automata to eventually synchronize their velocities resulting in a 
collective swirling motion of the entire cluster (Fig. 9b). In absence of confinement, 
our system of self-propelled ellipses shows the typical flocking behavior of the Vicsek 
model (Vicsek et al. 1995), so that for large density and small noise, the particles 
organizes in lanes or coherently moving subunits (see Fig. 10 and Supplementary 
Movie 6). We note that both collective swirling stasis originates from the interplay 
between self-propulsion, particle geometry and confinement and do not occur in 
systems without boundary. 

Our model described by Eq. (4.1) illustrates the origin of the three observed 
behaviors in a broader context. Analogously to the simple self-propelled particles, 
BBots tend to migrate to the boundary, which depending on the local density of the 
BBots and their angular velocity, can either play the role of an obstacle that causes 
the objects to jam, or a confining channel that collects and aligns the BBots into a 
coordinated moving cluster. For a cluster of walkers at the boundary, each BBot in 
the cluster is trapped by its neighbors and cannot escape. As their angular velocity 
increases, they can exert a sufficient torque on their neighbors to push them aside 
and escape from the cluster, consistent with the observation that the number of 
BBots required for jamming decreases with their angular velocity, as shown in the 
phase diagram in Fig. 6. Bbots that are spinners have a relatively small translational 
velocity and so are easily trapped by their neighbors at the boundary once their 
density is large enough and they unable to reverse direction and escape. However, 
the finite spinning torque leads to a global tilt of the BBots leading to a global 
swirling motion of the entire cluster along the edge of the arena. 
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Figure 9. Numerical simulations of collective behavior of Bbots. As in Fig. 5 the columns 
summarize the three behaviors observed in the experiments with BBots: (a) disordered 
(random) motion of spinners at low density; (b) swirling motion of spinners at high density, 
and (c) stasis of walkers at high density. (al-a3) Numerical solution of Eqs. (4.1) in the 
random phase, with uja/vo = 0.05; (al) Instantaneous position of 5 spinners with a = 5°, 
(a2) the mean velocity vi and the mean speed V2 of the Bbots, (a3) the trajectory of the 
center of mass of the Bbots in physical space, showing random motion. (bl-b3) Numerical 
solution of Eqs. (4.1) in the random phase, with uja/vo = 0.03; (bl) Instantaneous position 
of 15 spinners with a = 5°, (b2) vi and V2 showing a non-zero value, (b3) the trajectory of 
the center of mass of the Bbots in physical space, showing the signature of the coordinated 
swirling. (cl-c3) Numerical solution of Eqs. (4.1) in the random phase, with ct;a/vo = 0; (cl) 
Instantaneous position of 15 walkers with (c2) vi and V2 of the Bbots, (c3) the trajectory 
of the center of mass of the Bbots in physical space, showing no motion, i.e. stasis. For 
all simulations, we chose the ellipse aspect ratio to be 5, while the other parameters are 
kia/vQ — 10, foa^/vo = 1 and C'^/'^o = 27r. The equations are integrated via a four-step 
Runge-Kutta algorithm with time step At = 0.001. 

5. Discussion 

While the geometric structure of the clusters depends significantly on the shape 
of the particles, the occurrence of these three collective behaviors observed in the 
experiment is rather general. To demonstrate this we have run an additional set 
of simulations using self-propelled polar disks in place of elliptical particles (see 
Appendix B). The disks' dynamics is also dictated by Eq. (4.1), in which the 
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Figure 10. Example of flocking in a group of 100 ellipses from a numerical solution of Eqs. 
(4.1) on a periodic square domain of size Lja— 100. The right panel shows the evolution 
of the order parameter v\ and V2 as a function of time. The parameter values used are: 
cj = 0, = 1, /ci = 10, 0.2, C = 1- 



physical torque is Mij now replaced by a generic aligning interaction of the form: 
Mi^ = sm{Oj — 6i)/Ni (which has no real physical basis in our system). With this 
choice, Eq. (4.1b) becomes a short-range version of the Kuramoto model for phase- 
synchronization in chemical and biological oscillators (Acebron at 2005). Analo- 
gously to the ellipses, the dynamics of the polar disks is also characterized by three 
regimes: random motion at low densities, jamming at the boundary for large densi- 
ties when the angular velocity of the disks vanishes, and the formation of a compact 
cluster that circulates along the boundary for large densities when the disks have 
a finite angular velocity (see Appendix B). 

The similarity between our experimental system and the two models described 
above suggests that the coordinated circulation and jamming in a system of con- 
fined agents is generic. This form of collective behavior relies on simple but crucial 
features of the individual agents as well as the environment: the ability to translate 
and rotate, and the ability to interact with each other and with the environment, 
here including the boundary and the local topography. While the spatial structure 
of the clusters crucially depends on the shape and the packing properties of the 
particles, their collective motion is very robust, and depends on simple non-specific 
principles. 

In living systems, where similar behaviors such as the density-driven transitions 
are seen in confined Macrotermes michaelseni (Turner 2011), they have been linked 
with insect cognition and social interactions. Our study suggests that particle mo- 
tion, shape and spatial interactions are sufficient and might in fact play equivalent 
roles. In a biological setting such as termite swarms, one might test these ideas by 
controlling the confinement of termites by varying the substrate curvature and slip- 
periness, gluing circular discs on their backs to make the interactions more isotropic, 
etc. In an artificial setting, the collective abilities of spinner and walker BBots to 
convert environmental interactions into dynamical behavior may allow us to explore 
functional swarms that can search and sense environments. For example, they have 
the ability to sense substrate roughness by slowing down, and they can search to- 
pography (curvature) in massively parallel ways, using mechanical intelligence, and 
suggesting the use of these automata as fast, cheap, leaderless explorers. 
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Appendix A. Locomotion of an individual BBot 

The principle of motion of a single BBot, as inferred from high speed videos, relies 
on the sequence of events illustrated in Fig. 3. At each cycle of the eccentric motor 
the following sequence of events takes place: 1) the bristles bend as they are loaded 
by a force F = Mg + mruj'^^ where Mg is the weight of the BBot, m the eccentric 
mass, r the lever arm and uj the angular frequency of the motor. 2) While the 
eccentric mass rotates, the load on the bristles decreases; this causes the bristles 
to unbend. 3) The unbending bristles slip on the underlying substrate, producing 
a forward displacement of the object in the horizontal direction. 

A quantitative description of the gait reduces to calculating the horizontal dis- 
placement Ax of the bristles at each cycle of the eccentric motor. To accomplish this 
we ignore the collective dynamics of the bristles and focus on a planar description, 
replacing the rows of bristles as an ideal elastic rod subject to periodic tip-load 
acting in the y direction: 

Fy = W = Mg ^ mruj'^ sincjt . (A 1) 

When the eccentric mass is oriented with its axis of symmetry toward the negative 
^/-direction, the load is I^max = Mg-\-mruj^ and the bristles are maximally deflected. 
As the eccentric mass moves away from the vertical direction, the bristles start to 
recoil and their tip slides on the substrate. The sliding tip of the bristles is subject 
to a dynamic frictional force acting along the x-direction: 

F^ = fiW, (A2) 

where /i is kinetic friction coefficient. The unbending of the bristles terminates when 
the eccentric mass is oriented along the positive ^/-direction and the load is minimal: 
Wmin = Mg — mruj'^. 

To make progress we assume that the inertia of the bristles is negligible so that 
bristle deflection and sliding occurs quasistatically. This implies that the conforma- 
tion of the bristles is, at any time of the gait cycle, in equilibrium with the external 
load and the frictional force acting respectively on the y and x direction. Under 
this assumption, the shape of the bristles is governed by the classical equilibrium 
equations of an ideal elastic beam: 

Fs^K = 0, Ms^txF = 0. (A3) 

F and M are respectively the force and torque per unit length and the subindices 
denote a derivative with respect to the arc-length s of the beam. The tangent vector 
t of the bristles is given by: 



t = sinOx — cos y , 
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Figure 11. Schematic representation of BBot principle of motion. The two rows of bristles 
are modeled as a single elastic beam whose free end is subject to two forces: the time 
dependent weight W — Mg + mru)^ sin ut acting on the positive ^-direction and a kinetic 
frictional force acting on the positive x-direction. The gait cycle is assumed quasistatic so 
that the shape of the bristle is, at any time, in equilibrium with the applied forces. 

where is the angle formed by the bristles with the vertical direction. Finally, K 
is the external force acting on the tip of the bristles, thus: 

K = W{/^x^ y)6{s - L) , (A 5) 

where L is the length of the bristles and the delta function reflects the fact that 
the force is applied at the tip. Integrating the force equation and replacing it in the 
torque equation gives: 

Ms = WsinOz^ fiWcosOz. (A 6) 

The torque M acting in the beam and its curvature n = Og are related by the Euler- 
Bernoulli constitutive equation M /EI = —nh^ where h is the binormal vector of 
the beam (£ in this case) and EI is its bending rigidity (with E the Young modulus 
and / the area- moment of inertia). This yields a single differential equation for the 
angle 0: 

ElOss^WsmO^liWcosO = 0, (A 7) 

with boundary conditions: 

(9(0) = a, 6>,(L) = 0. (A 8) 

These are the typical boundary conditions of a cantilever beam, with one end fixed 
at an angle a and the other free of torques. It is convenient to work with dimen- 
sionless quantities, by rescaling the arc-length with the total length of the bristles: 
t = s/L, so that Eq. (A 7) can then be recast in the form 

^u + fe'sin(^ + (^) =0, (A 9) 



where: 



= j^j + 5 ^ = arctan ^ . 



(A 10) 
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Figure 12. Example of configurations of the deflected bristles obtained from Eqs. (A 20) 
for various ^Y/l"^ values. The orientation of the bristles at the supported end is a = 20° 
and the friction coefficient is /x = 0.1. 

Eq. (A 9) can be integrated exactly in terms of Jacobi elliptic functions to yield 

sin^(6> + (/9) = msn{k{t- 1) K,m) , (Alia) 

co^\{e^Lp) =dn(/c(t-l) + K,m) , (A lib) 



Ot = 2mk cn {k{t - 1) + i^, m) , 



(A 11c) 



where Eqs. (All) follow using standard techniques (see for instance Davis 1960) 
and we use the standard notation for elliptic functions and integral, i.e. given the 
elliptic integral of the first kind: 



u = F{(j)^ m) 



dt 



\/l — m? sin^ t 



(A 12) 



with < m? < 1 the elliptic modulus and (j) is the Jacobi amplitude: (j) = am(ii, m). 
From this it follows that 



sn(ii, m) = sin (/) , cn(ii, m) = cos , dn(ii, m) = \J 1 — m? sin^ (j) . (A 13) 

Finally, the quantity K in Eqs. (All), is the complete elliptic integral of the first 
kind: K = F(|^, m). This enforces the boundary condition at the free end t = 1: 

et{l) =2mkcn{K,m) =0. (A 14) 

The elliptic modulus m, on the other hand, is obtained from the boundary condition 
at the fixed end through Eq. (A 11a): 

sin ^{a -\- (p) = msn{K — k^m) . (A 15) 

With the solution (A 11) in hand, we can now construct a parametric equation for 
the shape of the deflected bristles by integrating the tangent vector t: 

r{s')= I dst{s)=f ds{smOx-cosOy). (A 16) 

Jo Jo 

In order to use Eq. (A 16) we first set A = sm^{6 -\- cp) and B = cos ^{6 -\- (p) and 
note that: 



cos = 2AB sin (p + {B'^ - A'^)cos(p, 



(A 17a) 



sin = 2AB cos cp - {B'^ - A^)sm(p, 



(A 17b) 
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Figure 13. The step size Ax/L as a function of the ratio M/m (left) and the bristle angle 
a from Eqs. (A 20). The parameters used for the plots are g/ruj'^ — 1 and 



= 0.1. 



and that the integrals of terms containing A and B are given, up to a constant, by: 



/ 



dtAB 



m 



cn {k{t- l)^K, m) 



I 



dt {B^ -A^) = -t^-E (am {k{t - 1) + i^, m) , m) 



where E is the elliptic integral of the second kind, defined as: 

E {(j), m) = f dt Vl-m'^ sin^ t . 
Jo 

Then, combining Eqs. (A 17) and (A 18) we obtain the following parametric expres- 
sion of the coordinate of the bristles: 



(A 18a) 
(A 18b) 

(A 19) 



x{t)/L = xq/L ■ 



2m 



cn (r, m) cos (p ■ 



t — - E (am (r, m) , m) 
k 



sin if , (A 20a) 



2777/ 

y{t)/L = yo/L^ cn (r, m) sin (/p ■ 



t — — E (am (r, m) , m) 
k 



cos (p , (A 20b) 



where we have called r = k{t — 1) -\- K for brevity. The integration constants xq 
and yo are set so that x(0) = y{0) = so that 

2 2m 
xq/L = — E (am (K — k^m) , m) sin ip H — — cn (K — /c, m) cos (f> , (A 21a) 
k k 

2 2m 

yo/L = - E (am (K — k^m) ^ m) cos ^ — cn (K — /c, m) sin , (A 21b) 

k k 

Eqs. (A20)-(A21), and (A 15) along with the definitions (A 10) give the shape of 
the bristles. In Fig. 12 we show a sequence of typical configurations obtained from 
this solution for various values of 



w 



Given the shape of the bristles, the step size Ax of a Bbot associated with 
each gait cycle is dictated by the position of the tip of the bristles. The latter 
can be obtained from Eqs. (A 20) by setting t = 1 and noting that cn(i^, m) = 0, 
am(iC,m) 7r/2. Then 



x(l) = Xq ^ Lsm.(p) 



2E{m) 



(A 22) 
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where E{m) = ^(|^, m) is the complete ehiptic integral of the second kind. Because 
of our quasistatic approximation, x(l) depends exclusively on the applied load and 
so the step size is simply given by the difference in the position of the tip associ- 
ated with the maximal (eccentric motor in the negative ^/-direction) and minimal 
(eccentric motor in the positive ^/-direction) load. Using the definition 

kl = ^^^^2 \/l + (A23) 
we can finally express the step size in the form: 



Ax = Aq — 2L sin(arctan /i) 



(A 24) 



where m± = m{k±) and Aq = xo(/c+) — xo(/c_). The model is valid only as long as 
/c?. > 0, which implies Mg > mguj'^. For light BBots, where this condition does not 
hold, locomotion is complicated by the fact that when the eccentric mass is oriented 
along the positive ^/-direction, there is an upward directed force that makes the BBot 
lose contact with the substrate. The resulting jumping motion then couples with 
the dynamics of the bristles making the gait cycle intractable with the methods 
used here. Fig. 13 shows a typical step size obtained from Eqs. (A 20) as a function 
of the ratio M/m and the bristles inclination angle a. 

Our analytical results allow us to capture the qualitative aspects of the motion 
of a single Bbot and its dependence on the magnitude and frequency of the eccentric 
driving motor, as well as the dependence on the mass of the Bbot, the orientation 
and length of the bristles, consistent with experimental observations. An alternative 
analysis of the locomotion of an individual BBot was carried out by DeSimone & 
Tatone (2012) using methods of geometric control theory (see Alouges et al. 2008). 



Appendix B. Collective behavior of self-propelled disks 

In order to gain insight into the origin and the generality of the behaviors observed 
in our experiments and numerical simulations of interacting Bristle-Bots (BBots), 
we also compared the results with those obtained from the numerical simulation 
of self-propelled disk-like particles that are isotropic. The particles have both a 
positional degree of freedom given by their center of mass Vi and an orientation 
rii = (cos6>^, sin^^) with the position and the angle 6i that evolve according to 
Eqs. (4.1), but, in contrast with the case of elliptical collisions where there is a 
physical torque that causes alignment, Mij is chosen to be: 

M,j=sm{Oj-0,)/N,. (Bl) 

With this choice, Eq. (4.1b), is a short-range version of the Kuramoto model for 
phase-synchronization (Acebron et al^ 2005) and can serve as a rather general model 
for aligning interactions among self-propelled particles, although it has no direct 
physical origin. 

For our simulations, we assume that the particles are confined to a circular 
domain of radius R centered at the origin. The interaction between the particles 
and the boundary that takes place is assumed to have an identical form to that 
used in the simulations of the elliptical particles. 
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Figure 14. Example of treadmilling in a group of 15 polar disks obtained from a numerical 
solution of Eqs. (4.1). The top frames shows the regression of a tracer disk (labeled in 
yellow) from the leading to the trailing end of the cluster. The tracer disk remains in 
proximity of the trailing end for about a loop and then starts its progression toward the 
front along internal side of the cluster. The parameter values used are: cj = 0.05, vq = 1, 
ki = 10, = 3, /C3 = 0.1, a = 1, i? = 11, ^ = 2a, C = 57r/2. 




400 600 
Time 



1000 



Figure 15. Example of breathing in a group of 300 polar disks obtained from a numerical 
solution of Eqs. (4.1). The right panel shows the evolution of the order parameter vi and 
V2 as a function of time. The dips in the V2 trajectory correspond to the configurations 
where the disks are densely packed at the center of the arena (pannel a in the left side 
of the figure), while the peaks in the vi trajectory denote the bursts as a consequence 
of which the disks suddenly migrate to the boundary. The parameter values used are: 
uj 0.1, vo = 1, /ci = 10, fe 0.2, k3 0.1, a = 1, = 28, ^ = a, C = vr. 



While the spatial structure and packing properties of the clusters depends on 
the details of the system, and in particular on the shape of the particles, the oc- 
currence of the coordinated behaviors (swarming, swirling and stasis) appears to 
be a very robust property of systems of self-propelled agents in a confined space. 
These behaviors are not sensitive to the presence of inertia [this is present in the 
experiment and is neglected in Eqs. (4.1)] or to the shape of the particles and the 
precise form of the aligning torque Mij. However, there are features that do depend 
on the details; a most interesting example is what we term treadmilling, observed 
in the numerical simulation of both ellipses and disks. Fig. 14 shows an example of 
treadmilling in a group of polar disks, wherein particles move through the cluster 
and eventually leave it, only to join it later at the other end. For very large den- 
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sities, the self-propelled disks also exhibit a breathing mode^ in which the particles 
periodically migrate from the center to the boundary and vice-versa by mean of 
sudden bursts (Fig. 15) reminiscent of those observed in excitable active systems 
(see Giomi et al 2011 and 2012). 
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